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COPYRIGHT (c) 1978, 1980, 1982, 1984 BY 
DIGITAL EQUIPMENT CORPORATION, MAYNARD, MASSACHUSETTS. 
ALL RIGHTS RESERVED. 


THIS SOFTWARE IS FURNISHED UNDER A LICENSE AND MAY BE USED AN 
ONLY IN ACCORDANCE 

INCLUSION OF THE ABOVE COPYRIGHT NOTICE. THIS SOFTWARE OR ANY OTHER 
COPIES THEREOF MAY NOT BE PROVIDED OR OTHERWISE MADE AVAILABLE TO ANY 
OTHER PERSON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY 
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D COPIED 
WITH THE TERMS OF SUCH LICENSE AND WITH THE 


THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE 
epokat ite NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT 
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FACILITY: Language support library - user callable 
ABSTRACT: 


D COMPLEX*16 base to INTEGER*4 power. 
Floating overflow can occur. 
Undefined exponentiation can occur if 
base = (0.,0.) and exp <=0 


VERSION: 1 
HISTORY: 
AUTHOR: 
Steven B. Lionel, 27-July-1979 


J 9 
T CDJ = D COMPLEX*16 ** INTEGER*®4 power routin 16-SEP-1984 01:55:4 AX/VMS Macro v04-00 Pa 
9t35 HISTORY ; Detailed Current Edit History mitt 1 91:33:28 YRTARTL Se OTSPOWCDJ .MAR; 1 - 3) 


2 -SBTTL HISTORY ; Detailed Current Edit History 
; 8 = Adapted from OTSS$POWCJ version 1-003. SBL 27-July-1979 
1- - correct bug in ap Ae tor yndet ines result with negative powers. 
SPR 11-35262 SBL 22-Jan-198 
1-003 = Use general mode addressing. SBL 30-Nov-1981 
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.SBTTL DECLARATIONS 
$ INCLUDE FILES: 


; EXTERNAL SYMBOLS: 


tine 
on routine 


: MACROS: 


: EQUATED SYMBOLS: 
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OWN STORAGE: 
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' : PSECT DECLARATIONS: 
Be -PSECT _OTSSCODE PIC,SHR,LONG,EXE ,NOWRT 


Hi 3 program section for OTS$ code 
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cae SBTTL OTSSPOWCDJ_R3 = D COMPLEX#16 #* INTEGERS4 
; FUNCTIONAL DESCRIPTION: 


D COMPLEX*16 result = D COMPLEX®*16 base ** INTEGER*4 exponent 
The COMPLEX result is given by: 
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base exponent resuct 

any >0 PRODUCT (base * 2**i) where 
i is each non-zero bit in 
exponent. 

(0., 0.) <=0 Undefined exponentiation. 

not (0., 0.) <0 PRODUCT (base * 2**i) where 
i is each non-zero bit in 
‘exponent:. 

not (0., 0.) =0 (1.0, 0.0) 


Floating overflow can occur. 

Undefined goes occurs if base is 0 and 

exponent is 0 or negative. 
CALLING SEQUENCE: 

result.wdc.v = OTSSPOWCDJ_R3 (base.rdc.v, exponent.rl.v) 
INPUT PARAMETERS: 

base = 4 ; D COMPLEX*16 base passed ¥ VALUE ! 

exponent = 20 ; Longword integer exponent by value. 


IMPLICIT INPUTS: 
NONE 


00000004 
00000014 


OUTPUT PARAMETERS: 
NONE 
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IMPLICIT OUTPUTS: 
NONE 


FUNCTION VALUE: 


THE D COMPLEX*16 result is returned in registers RO-R3. 
This is a violation of the VAX calling standard, but is 
excused for compiled code support routines. 


SIDE EFFECTS: 


Modifies registers RO-R3! 

SS$_FLTOVF = tanice’t gverficw : 
SIGRALs MTHS_UNDEXP (82 = * UNDEFINED EXPONENTATION') if 
base is 0 and exponent is 0 or negative. 
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ASHQ 
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sENTRY OTSSPOWCDJ_R3, “M<R4,R5, R6, of? »R8> 


MOVQ base (AP) 
MOva peeesBiAb) 
MOVL exponent (Ab), RB 


epee’ #0, RB. EVEN 


MOVD 

MOVD Re. 4 

ROTL #1, RB, RB 
BEQL oné 


Dd 
BRB SQUAR1 


MOVD #1, RO 


ROTL #-1, R8, RB 
BNEQ SQUARI 


TSTD 
BNEQ DONE 
TSTD R6 
BNEQ DONE 
#15, #1, RO 
ASHQ #1 R2 
MOVZBL aMTHsK K_UNDEXP, -(SP) 
CALLS , G*ATHS$SIGNAL 
RET 


ASHL #-1, R8, R8 
= square current base 


“ MULD3 -R4, R6, =(SP) 


MULD2 R4, RG 

MULD RG, R6 

SUBD R4 

ADDD Ree, ® (SP)+, R6 
BLBC RB, Sau 
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sable integer overflow 


> ROR gets COMPLEX base 
; RB = Longword exponent 


8 = | exponent } 


; pron h if even and clear low bit 


RO-RS = initial result 
R8 = unsigned_exponent / 2 


; done if exponent was 1 
; else use rest of exponent 


RO-R3 = initial result 

(1.0, 0.0) 

R8 = unsigned_ exponent / 2 
branch if exponent not 0 
exponent was 0, text RP(base) 
done if non-0, answer is 1.0 


: IP (base) better nok be zero 


t isn't return 1 


return RO-R3 = reserved operands 
FORTRAN ~_— ?<%} 


ert to 32-bit condition code 
and SIGNAL OMTHS UNDE XP 


R8 = ‘reduced exponent: / 2 


; +t By tmp = RP(base)*IP(base) 


R4-R5 = RP(base)** 
R6-R7 = IP(base)** 
R4-R5 . RP(base)**2 - IP(base)**2 


2* (RP (base) *IP (base) ) 


; branch if next exponent bit is 0 


= partial result * current power of base 


MULD RO, R6, -(SP) 


MULD R 
MULD3 = R2, RG. 5 (SP) 
SUBD2 (SP)+ 

MULD R 

ADDD2 (SP)+ 

ASHL #-1 gg? RB 


BNEG _ SQUARI 
TSTL exponent (AP) 


RO-R1 = part) * RP 
(SP) = mi = 1P(part) iBtbas 

RO-R1 = R c,iPtoact) IP [P(part)*1P(base 
4 “R35 = IP(part)*RP (base) 
R2-R5 = IPC part) *RP(base) +RP(part)*1P (base 
R8 = ‘reduced exponent: 


loop if more exponent bits left | 
test exponent sign 
| 
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POWCDJ 
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RECIP 
UNDEF INED 
R2, -(SP) 
RO, =(SP) 

-(§P) 

#1, =(SP) 

#8, G*OTSSDIVCD_R3 


91:33:68 yaaivas Macro _V 


04-00 | 
4 MTHRTL.SRCJOTSPOWCDJ.MAR; 1 
done if positive | 
test RP(result) 
if non-0, OK to take reciprocal 
RP(result) was 0. test IP(result) 
undefined (0.0+0.01) ** <n 


; second arg pair is divisor 
; push (1.0,0.0) on stack 

; RO-R3 = reciprocal 

; result in RO-R3 
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OTSSPOWCDJ 255:4 AX/VMS Macro V04-00 Pa 7 
Symbol table 66-SEP-1984 atk 1:27: 4 MTHRTL.SRC €10 OTSPOWCDJ .MAR; 1 “t (5) 
BASE = 00 
DONE ¢ R 4 
EVEN R 1 
EXPONENT = ‘ 
MTHSS$SIGNAL eetereee = §6X 0 
MTHSK_UNDEXP ereeeeee YX 0 
OTSSDIVCD_R3S ereeeeer = §6X 0 
OTSSPOWCDJ_R3 O90000A1 00 ne 01 
POWCDJ 4 
RECIP 44 ar R 1 
SQUAR 0000 rt R 1 
SQUAR1 0000 3 R 1 
UNDEF INED OOGOO0SA R 1 

eeeceeoeesooeoeoes + 

: : Psect synopsis ; | 
PSECT name Allocation PSECT No. Attributes 
»— OS , 00000000 <( 0.) 00 ¢ 0.) NOPIC USR CON ABS LCL ee NOEXE NORD NOWRT NOVEC BYTE 
_OTSSCODE OOOOO0A2 ( 162.) O1¢ 1.) PIC USR CON REL LCL SHR EXE RD NOWRT NOVEC LONG 

tec eee wesc e eee eee eee cee ces 
i Performance indicators ! 

Phase Page faults CPU Time Elapsed Time 


00:00:00.97 | 
| 


Initialization 9 00:00:00.09 

Command processing 124 00:00:00.45 00:00:02.59 
Pass 76 00:00:00.60 00:00:01.93 
Symbol table sort q 4 SB 00:00:00.04 
Pass 2 5 00:00:00.4 00:00:01.53 
Symbol table output 00:00:00.02 O33 +20 «Bs 
Psect synopsis output 00:00:00.01 00:00:00.0 

Cross-reference output OF RD py oe 00:00:00.00 
Assembler run totals 28 00:00:01.65 00:00:07.10 


The working rs Limit was 900 pages. 
3177 bytes (7 pages) of virtual memory were used to buffer the intermediate code. 
There were 10 pages of symbol table space allocated to hold 13 non-local and 1 local symbols. 
19 source Lines were read in Pass 1, producing 11 object records in Pass 2. 
pages of virtual memory were used to define 0 macros. 


See tener een ee corona ce see a} 


Macro library name Macros defined | 
“S255SDUAZB:(SYSLIBISTARLET.MLB;20 0 SSS 0 | 
0 GETS were required to define 0 macros. 

| 


There were no errors, warnings or information messages. 
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MACRO/ENABLE=SUPPRESSION/DISABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:OTSPOWCDJ/OBJ=OBJ$:OTSPOWCDJ MSRC$:OTSPOWCDJ/UPDATE=(ENHS$:OTSPOWCDU) 


oTsS CDJ 
VAX=-11 Macro Run Statistics 


244 AH BT13A-SE : ENT CORPORATION 
VAX/VMS V4.0 ND PROPRIETARY 


